Model evidence of a superconducting state with a full energy gap in small cuprate 

islands 
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We investigate subdominant order parameters stabilizing at low temperatures in nano-scale high- 
T c cuprate islands, motivated by the recent observation of a fully gapped state in nanosized 
YBa2Cu307-5 [D. Gustafsson et a/, Nature Nanotech. 8, 25 (2013)]. Using complementary quasi- 
classical and tight-binding Bogoliubov-de Gennes methods, we show on distinctly different properties 
dependent on the symmetry being d x i_ y i + is or d x i_ y i + id xy . We find that a surface-induced 
d x 2_ y 2 + is phase creates a global spectroscopic gap which increases with applied magnetic field, 
consistent with experimental observation. 

PACS numbers: 74.20.Rp, 74.50.+r, 74.72.Bk 



It is well established that high-temperature cuprate su- 
perconductors have a dominantly d x 2_ y 2-wdive order pa- 
rameter symmetry [l|, |2| , but the existence of a subdomi- 
nant symmetry has also long been considered. Of special 
interest are order parameters such as d x 2_ y 2 -\-id xy (d\ + 
id^) and d\ + is, which fully gap the Fermi surface and 
break time-reversal (T) symmetry [3|— ls[ - Experimental 
data have been contradictory, invoking large imaginary 
subdominant orders to explain tunneling experiments in 
YBa 2 Cu 3 7 _ 5 (YBCO) 0,0 and La^Sr^CuC^ @ or 
thermal conductivity in Bi 2 Sr 2 CaCu208 [lQj while, on 
the other hand, only very small imaginary components 
seem compatible with the absence of any measured spon- 
taneous magnetization [2|, lllMl4 1 . 

The possibility to find a T-symmetry breaking state 
is enhanced if the dominant di-wave order parameter is 
locally reduced by e.g. surface scattering. A nano-scale 
island of a curate superconductor is, by virtue of its large 
surface-to-area ratio, thus an ideal candidate to search 
for this elusive state. Very recently, an even/odd par- 
ity effect was reported in YBCO single-electron transis- 
tors (SETs), signalling a fully gapped low-temperature 
superconducting phase in nano-scale YBCO [15]. For 
a pure di-wave superconductor there are always low- 
energy quasiparticle states available at the nodal points, 
into which the added charge in a SET can relax. Un- 
less the nodal quasiparticles states are lifted by a spec- 
troscopic energy gap E g , a parity effect should not be 
present in YBCO-SETs. The experimental data show 
Eg « 20 — 40 /ieV, making E g three orders of magnitude 
smaller than the gap A^ ~ 20 meV in YBCO |15J . 

The YBCO SETs studied in Ref. [15] have sizes of order 
0.2 x 0.2 x 0.1 /im 3 , making the energy- level spacing S s 
of the nodal quasiparticles a candidate for the observed 
gap Eg. While 5 S ~ 10 /ieV for an infinite mean- free 
path, surface disorder gives a reduced S s ~ 0.01 /xeV 



<C Eg [16j. However, energy level spacing due to finite 
size is incompatible with the finite onset voltages of the 
SET current measured over a full gate charge period [17j . 
The measured E g is thus of superconducting origin and 
gap the whole Fermi surface. Given the parent di-wave 
order parameter of cuprates, finding E g demonstrates the 
presence of a complex order parameter, which gaps the 
nodal di-wave spectrum. 

The aim of this Letter is to establish which fully 
gapped superconducting state nucleates at low tempera- 
tures in nano-scale cuprate islands. Specifically, we focus 
on three particularly illuminating experimental results in 
Ref. [15]: a) The spectral gap E g « 20 - 40 /ieV, b) the 
energy level spacing 5 S ~ 1 /ieV above E g , and c) E g 
increases with applied magnetic field in the range 0-3 T. 
We use complementary quasi-classical and tight-binding 
Bogoliubov-de Gennes methods to show that the d\ + is 
and d\ + id2 states have distinctly different properties in 
small cuprate islands. While a subdominant s-wave com- 
ponent appears at the surface due to disorder suppression 
of the di-wave, a <i 2 component can only nucleate in the 
interior of the island. Furthermore, the d\ + is state has 
a finite gap set by the decay of the s- wave order into the 
center of the island, whereas the low-energy spectrum of 
the d\ J rid2 state is determined by finite size quantization 
of its chiral surface states. Thus, the d\ + is state has 
a low-energy spectrum similar to that of a conventional 
s- wave superconductor, qualitatively satisfying points a) 
and b), while the d\ + zd 2 state has equally spaced low- 
energy levels. We also find the magnetic field dependence 
to only be consistent with the d\ + is state. 

The necessary attraction in a subdominant pairing 
channel for a gapped state is an inherent feature of boson- 
mediated pairing, such as a spin-fluctuation model rele- 
vant for high-T c cuprates [18I422J. The relative strength 
of the attraction in different paring-symmetry channels 
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FIG. 1: (Color online) Subdominant s-wave order, (a): QC s-wave component A s in the center of the island (o) and nodal energy 
gap Eg (x) of a pair-breaking [110] surface with disorder as function of slab length L/^o for T c , s /T c = 0.001 (magenta), 0.01 
(blue), 0.1 (green) (increasing values) at temperature T = 0.01T C . Dashed lines are A s (center) /Ad 1 = 4 A s (surface) /(L/£o), 
dotted lines are E g /Ad 1 = 48 A s (surface) /[(l — log(T C)S /T c ))(L/£o) 2 ]. (b): BdG disorder averaged energy gap E g as function of 
A s in the center of the grain (solid) and at the surface (dashed) for T c , s = 0.5T C (black, o), 0.25T C (red, x), and 0.1T C (green, 
D). Dotted line is E g = 2 A. (c): BdG disorder averaged LDoS in the center of the grain (solid) and at the surface (dashed). 
(d): BdG eigenstate spatial density for one disorder configuration with L = 40 (black dots marking removed sites) averaged 
over the 4 lowest energy states (left) and disorder averaged (40 configurations) for the lowest energy state (right). White = 
zero, black = 0.01 (left) or 0.005 (right) states per unit cell, (e): BdG disorder averaged level spacing S s — E n +i — E n ratio 
to first energy value E\ for n — 1 (solid) and n — 2 (dashed) as function of grain size L. (f): Evolution of the QC energy gap 
Eg (black, x) and A s on the surface (cyan, □) with magnetic field, Bq = <J>o/vr£oAo ~ 2.2 T (derived assuming £o = 2 nm, 
Ao = 150 nm for YBCO) for a L = 40 £o slab with surface disorder modeled by a thin layer (~ 0.2£o) with a graded impurity 
concentration [25|] and T C)S = 0.065T C . 



is sensitively dependent on the shape of the pairing sus- 
ceptibility and the band structure in vicinity of the Fermi 
level. In small cuprate islands, with spatial dimensions 
down to 50-100 times the superconducting coherence 
length £o i finite size effects might become important such 
that the pairing interactions deviate from the bulk val- 
ues. Here we will not further elaborate on this but in- 
stead focus on the consequences of an assumed subdomi- 
nant pairing by following Refs. [3|,|23j to generate generic 
phase diagrams for the d\ + is and d\ + id2 states. We 
quantify the strength of the subdominant cfe- or s-wave 
pairing simply by their bare bulk transition temperature 
T c d 2 / s , measured relative to T c , the transition tempera- 
ture into the d\ state. 

To treat the appearance of subdominant orders in a 
cuprate island we employ two complementary methods. 
First we calculate the nucleation of subdominant order 
parameters in [100] and [110] surface slabs within the 
weak-coupling quasiclassical (QC) approximation for su- 
perconductivity (see e.g. Ref. [24]). The advantage of 
this approximation is that we can self-consistently com- 
pute the order parameter mean-fields in restricted ge- 
ometries of realistic size, including effects of both single- 
impurity scattering, impurity self-energies, and surface 
scattering. A self-consistent evaluation of both ther- 
modynamical properties, such as order parameter fields 
and magnetization, and the spectral properties is done 



as function of subdominant pairing strength and applied 
magnetic field. Since the QC approximation is a leading- 
order theory in quantities such as l/(kp^o) ~ 0.1 — 0.2 for 
high-T c cuprates, it cannot, however, resolve effects of a 
finite level spacing S s . To complement these results, we 
therefore also study two-dimensional (2D) square lattice 
islands within the Bogoliubov-de Gennes (BdG) frame- 
work [26]: 

#BdG = J2 ^j c I c j + J2 A di(hj)[4t c ^ - 4i c U (*) 

M> (ij) 

+ X>*w4r c l + E A d 2 {hj)\4Ai - c l c U + Hc - 

Here we use a band structure parameterization t^ rel- 
evant for a 2D model of high-T c cuprates [27|, includ- 
ing up to next-next nearest-neighbor hopping. The 
dominant d\ -pairing is implemented on nearest-neighbor 
bonds [28[ with the (mean-field) order parameter calcu- 
lated self-consistently using A^ = — Vd 1 (c^c J -^ — c^Cji}, 
with Ad 1 (x) and Ad x (y) treated independently. Due to 
computational demands limiting our grain sizes, we set 
Vd x = 0.455 eV giving coherence peaks at energies 5 
times larger than in YBCO @. Thus, while the BdG 
approach accurately captures the short £o in cuprates, 
we are instead limited by studying smaller grains and 
stronger superconductivity than found experimentally. 



Subdominant s- wave pairing is introduced through a neg- 
ative on-site potential V 3 , with self-consistency equation 
A s = —V 3 (ci±Ci^). For subdominant d^ pairing we use 
a next-nearest neighbor pair potential Vd 2 analogous to 
Vd-L but with imposed ^2 symmetry. We study islands 
with an overall [100] square shape with sides as large 
as L = 50 unit cells, thus L/^o is comparable to recent 
experiments [15J. We introduce surface disorder by ran- 
domly removing up to 25% of the surface atoms in the 
three outermost surface layers, see Fig. QJd). While no 
disorder average errorbars are displayed in the data in 
Figs. 1-2, they are negligible for all but the smallest is- 
lands. We find no significant difference in energy levels 
and subdominant orders between these disorder configu- 
rations and those of a diamond shaped [110] island with 
moderate surface disorder, thus establishing the generic 
nature of our islands. 

Subdominant s-wave order. — Surface scattering 
severely suppress <i-wave order at generic surfaces. Using 
both the QC and BdG methods we find that, given any 
finite pairing interaction in the s-wave channel, an s- 
wave order parameter with relative phase tt/2 nucleates 
in the d\ -voided disordered surface region, producing a 
T-symmetry breaking d + is state. A s is constant in the 
surface region, only determined by the s-wave pairing 
strength. The surface s-wave state leaks into the center 
of the island, with A s (center) oc A s (surface) / L 2 for the 
experimentally relevant grain sizes. We have analytically 
verified this power law decay by expanding the free 
energy of the bulk superconductor in A s to second 
order while keeping the dominating <i-wave component 
Ad ± exactly. The correction to the free energy takes 
the form SF S = J ^x s fe, q y )A* 3 (-q)A 8 (q), where 
the kernel x(Qx,q y ) oc 1 + £ s \q x - q y \ + ^ s \q x + q y \ 
reflects the strong anisotropy of the quasiparticle 
dynamics, and £ s ~ V s /Ad 1 is the s— wave coherence 
length. Accordingly, we find that, in agreement with 
the numerics, far away from a long boundary located 
at x = 0, the s-wave component should decay as 

Both the QC and BdG results have a finite gap E g 
in the energy spectrum, due to the finite s-wave order 
gapping the d\ nodes. The QC spectral gap E g decays 
as A s (surface)/L, as seen in Fig. QJa). This scaling is 
a consequence of the Fermi-surface position (pf ) depen- 
dence of the effective coherence length, being very long 
in the nodal direction of a <i-wave superconductor. In 
Fig. 1(b) we see that the BdG E g for small islands de- 
pends monotonically on A s (center), largely independent 
on both island size and s-wave pair potential. How- 
ever, the value of E g is much larger than the local result 
Eg = 2 A s (center), because of the large A s (surface), and 
for larger islands the surface order parameter ultimately 
determines E g . In Fig. 1(c) we plot the BdG local den- 
sity of states (LDoS) in both the center and surface re- 



gions which shows that the gap is a global property of 
the island, only the LDoS above the gap is position de- 
pendent. The non- localized spatial density of the lowest 
energy eigenstates in Fig. 1(d) further cements the fact 
that the gap depends globally on the subdominant order 
parameter. Further focusing on the lowest energy levels 
and their spacings, we plot in Fig. 1(e) the ratio of 5 S to 
the lowest energy level E\ = E g /2. This ratio decreases 
with increasing T c ^ s for two reasons: first, E g increases 
and second, S s decreases due to a more pronounced coher- 
ence peak [see Fig. 1(c)] above the s-wave nodal gap. The 
low-energy spectrum of a cuprate nano-scale island with 
a subdominant s-wave state nucleated only on the surface 
is thus essentially that of a conventional s-wave supercon- 
ductor. We find that in order to achieve S s /Ei ~ 0.05, 
as found experimentally |15j |. T C)S ~ 0.5T C for Ad t both 
2.5 and 5 times larger than the experimental value. This 
indicates that the s-wave pairing might be strong close 
to the surface of nano-scale cuprate islands, although the 
energy gap is still very small, since T C}S ~ 0.5T C only 
supports a surface s-wave state. Finally, in Fig. 1(f) we 
plot the QC magnetic field dependence for a surface dis- 
ordered slab. A s (surface) grows with magnetic field and 
we see that E g closely tracks this behavior, consistent 
with experimental results [15j. 

Subdominant d xy -wave order. — While s-wave order 
survives significant disorder, that is not the case for ^2 
subdominant pairing. We find using both the QC and 
BdG methods that a ^2 component nucleates only away 
from the surface, see Fig. 2(a), and then with a tt/2 rel- 
ative phase shift. We find that T Cj d 2 ^ 0.4T C is needed 
for a d\ + id2 state [23j , which requires an enhancement 
of the pairing attraction in the ^2-wave channel due to 
finite size effects. We will here not further discuss the 
probability of small dimensions increasing T c ^ 2 , but in- 
stead assume this can be the case and focus on the con- 
sequences. Any d\ + idi state, even with a very small 
c?2 component, hosts two chiral edge states [J, ||J |30|, in- 
dependent on the surface morphology. Fig. 2(b) shows 
how the surface states produce a constant LDoS in the 
surface region, due to their ID Dirac spectrum, and how 
this density leaks into the middle of the island in the QC 
results. The finite density of surface states in the cen- 
ter produces a very small hybridization gap as evident in 
the DoS very close to zero energy. We estimate the level 
spacing in these chiral modes to be S s ~ hv g Sk which 
is on the order of 10 fieV for an island circumference of 
500 nm. Here the gap- velocity v g = 2Ad 1 /hkY measures 
how Ad t opens in the node. In Fig. 2(c) we complement 
the QC result by plotting the lowest energy levels in the 
BdG results. Notably, they scale as L _1 ~ 5k, are inde- 
pendent of T Cj d 2 , and are equally spaced, i.e. S s « E\. 
These results all confirm that the lowest energy states in 
a d\ + id2 superconducting island are those of the chiral 
edge states. This is also evident from the spatial density 
of the lowest energy states shown in Fig. 2(d). The only 
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FIG. 2: (Color online) Subdominant cfo-wave order, (a): QC di-wave (dashed), and cfo-wave (solid) components across a 
L = 40^o [100] slab for T c?d2 = 0.5T C (red, x) and 0.47T C (cyan, o). (b): QC LDoS in units of normal state LDoS in the center 
(black) and at the surface (red). Inset shows a zoom-in at low energies, (c): BdG energy eigenstates times grain size L as 
function of L for E x (solid) and E 2 (dashed) for T c , d2 = 0.6T C (black, o), 0.5T C (red, x), 0.4T C (green, □), and 0T C (blue, o). 
(d): Same as Fig. 1(e) but for cfo-wave order, (e-f): QC band structure at the surface for B = (e) and B = 2Bo (f) with color 
scale showing the angle resolved LDoS, i.e <p P{ measures the position on the Fermi circle measured from the k x -axis. Dashed 
lines mark the d\ node. 



possibility of avoiding measuring the chiral edge states 
would be in small islands with 5 S > E g , where E g is the 
nodal gap in the center of the island. However, we see in 
Fig. 2(c) that even with T c ^ 2 = 0, the spectrum is still 
equally spaced, now from a finite size quantization in the 
d\ nodes. Both low-energy level spacings and eigenstate 
densities are thus distinctly different for a ^2 subdom- 
inant order compared to an s-wave order. In Fig. 2(e, 
f) we examine the magnetic field dependence of the QC 
band structure. The main effect of the magnetic field is 
to shift the zero energy momentum of the edge modes. 
Thus, 5 S will not change with field. We also note that the 
QC magnitude of the subdominant &<i component does 
not noticeably grow with magnetic field for fixed pairing 
strength. Thus, neither energy level spacings nor mag- 
netic field dependence of the d\ + id2 state is seemingly 
consistent with current experimental data |15j. We fi- 
nally note that self-consistent BdG results for both finite 
s-and ^2-wave pairing strengths, in general results in only 
one emergent subdominant order, despite their spatially 
separated nucleation regions, thus producing spectra sim- 
ilar to those already discussed. 

Spontaneous currents. — Both the d\ + is and d\ + id2 
state generate spontaneous currents due to T-symmetry 
breaking. The current experimental status regarding a 
superconducting state breaking T-symmetry taken from 
scanning- SQUID experiments [13| and /3-detected nu- 
clear magnetic resonance 14] put a strict upper limit 
on a subdominant order parameter A su b < 0.02A C ^ 1 [13j. 
Also the spontaneous magnetization is limited by an up- 
per bound of 0.2 G [14|]. Within the BdG framework we 
can calculate the quasiparticle currents by combining the 



charge continuity equation with the Heisenberg equation 
for the particle number (see e.g. Ref. [31]). In Fig. 3 
we plot the disorder averaged clock- wise surface currents 
along each four sides of an island. The d\ + id2 solution 
has a circulating surface current. While spontaneous, 
the current is not quantized, in agreement with d\ + id2 
superconducting graphene [32]. For d\ + is the surface 
current instead closes in small separate loops that form a 
staggered pattern of clockwise and anticlockwise current 
flow, with directions displayed in the inset. These local- 
ized current vortices will significantly reduce the mag- 
netic field associated with the spontaneous current, con- 
sistent with the present experimental situation [13|, LL4[ . 
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FIG. 3: (Color online) Disorder averaged surface currents in 
the clockwise direction summed over the 10 first surface lay- 
ers on each four sides (left, up, right, down) for L — 40 and 
T c , d2 = 0.5T C (thick red) and T c , s = 0.5T C (black). Verti- 
cal dashed lines mark the corners. Insets show schematically 
the surface current orientation. Black arrows indicate cur- 
rents present after disorder averaging, blue arrows indicate 
outermost surface currents prone to cancelation in disorder 
averaged results. 



In summary we have shown that the time-reversal sym- 
metry breaking d\ + is state in a cuprate island has a fi- 
nite energy gap, above which the subsequent energy level 
spacings are dense. This is very distinct from the d\ -\-id2 
state where the low-energy spectrum is determined by fi- 
nite size quantization of the two chiral surface states cir- 
cling the island. In an applied magnetic field, the energy 
gap increases sub-linearly in the d\ + is state, whereas 
there are no significant change in energy levels for the 
d\ + id2 state. Comparing to recent experimental results 
[15| our results indicate that a d\ + is state might be 
present in small YBCO islands. 
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